ALGEBRO-GEOMETRIC CONSTRAINTS ON SOLITONS WITH 
RESPECT TO QUASI-PERIODIC BACKGROUNDS 
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Abstract. We investigate the algebraic conditions that have to be satisfied 
by the scattering data of short-range perturbations of quasi-periodic finite-gap 
Jacobi operators in order to allow solvability of the inverse scattering problem. 
Our main result provides a Poisson-Jensen-type formula for the transmission 
coefficient in terms of Abelian integrals on the underlying hyperelliptic Rie- 
mann surface and an explicit condition for its single- valuedness. fn addition, 
we establish trace formulas which relate the scattering data to the conserved 
quantities in this case. 



1. Introduction 

Solitons are a key feature of completely integrable wave equations and there are 
usually two ways of constructing the iV-soliton solution with to respect to a given 
background solution. Both are based on fact that the underlying Lax operator is 
reflectionless with respect to the background, but has N additional eigenvalues. 
One is via the inverse scattering transform by choosing an arbitrary number of 
eigenvalues (plus corresponding norming constants) and setting the reflection coef- 
ficient equal to zero. The other is by inserting the eigenvalues using commutation 
methods. This works fine in case of a constant background solution and the eigen- 
values can be chosen arbitrarily. However, in case of a (quasi-)periodic background 
solution it turns out that the eigenvalues need to satisfy certain restrictions. This 
was probably first observed in [14] , where it was proven that adding one eigenvalue 
to the two-gap Weierstrass solution of the Korteweg-de Vries (KdV) equation pre- 
serves the asymptotics on one side, but gives a phase shift on the other side. The 
general case was solved in |S]. In particular, this shows that the eigenvalues and 
reflection coefficients can no longer be prescribed independently if one wants to stay 
in the class of short-range perturbations of a given quasi-periodic background. It 
turns out that these constraints are related to the fact that the resolvent set of the 
background operator is not simply connected in the (quasi-)periodic case. In this 
case we have to reconstruct the transmission coefficient from its boundary values 
on this non simply connected domain which is only possible in terms of multivalued 
functions in general, see [25j . Hence one needs to impose algebraic constraints on 
the scattering data to obtain a single-valued transmission coefficient. It seems that 
this was first emphasized in [4]. 

The aim of the present paper is to make this reconstruction explicit in terms of 
Abelian integrals on the underlying hyperelliptic Riemann surface for the case of 
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Jacobi operators (respectively the Toda equation). However, similar results apply to 
one dimensional Schrodinger operators (respectively the KdV equation). This will 
then allow us to derive an explicit condition for single-valuedness and to establish 
trace formulas which relate the scattering data to conserved quantities for the Toda 
hierarchy. In particular, these trace formulas are extensions of well-known sum rules 
(see e.g. [3], [IT], [15], [IT], [18], [26]) which have attracted an enormous amount of 
interest recently. 

To achieve this aim we will first compute the Green function, harmonic measure, 
and Blaschke factors for our domain. This case seems to be hard to find in the 
literature; the only example we could find is the elliptic case in the book by Akhiczcr 
P]. See however also [21], [22], where similar questions are investigated. 

2. Notation 

To set the stage, let M be the Riemann surface associated with the function 
^2g+2 0)) where 

2g+l 

(2.1) R2 9 +2(z) = [| [z Ej), E < E 1 < ■ ■ ■ < Eag+i, 

3=0 



g £ N. M is a compact, hyperelliptic Riemann surface of genus g. We will choose 
R 2 g^ 2 (z) as the fixed branch 

2g+i 

(2.2) <+ 2 W = - II V^E~ 

3=0 

where yC is the standard root with branch cut along (— oo,0). 

1 /2 

A point on M is denoted by p = (z, ±R 2g+2 (z)) = (z, ±), z S C. The two points 
at infinity are denoted by p = oo±. We use 7r(p) = z for the projection onto the 
extended complex plane C U {oo}. The points {(Ej, 0), < j < 2g + 1} C M are 
called branch points and the sets 

9 

(2.3) n±={(z,±< 2 (z))|zeC\E}cM, E=\J[E 2j ,E 2j+1 ], 

3=0 

are called upper and lower sheet, respectively. Note that the boundary of II± 
consists of two copies of S corresponding to the two limits from the upper and 
lower half plane. 

Let {cij,bj} 9 j =1 be loops on the Riemann surface M representing the canonical 
generators of the fundamental group ~k\ (M) . We require a,j to surround the points 
E 2 j-i, E 2 j (thereby changing sheets twice) and bj to surround Eq, E23-1 counter- 
clockwise on the upper sheet, with pairwise intersection indices given by 

(2.4) cij-oafc =bjob k = 0, ajob k = S jk , 1 < j, k < g. 

The corresponding canonical basis {Cj'}f=i f° r the space of holomorphic differentials 
can be constructed by 

■k 3 dn 



(2.5) C = EsCi) 



pl/2 
3 = 1 rt 2 ff +2 
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where the constants c(.) are given by 
The differentials fulfill 

(2-6) / Cfc = <5j,fc, / (k=T j} k, T j,k=Tk,j, l<j,k<g. 

J aj J bj 

For further information we refer to [6], [201 App. A]. 

3. ALGEBRO-GEOMETRIC CONSTRAINTS 

We are motivated by scattering theory for the pair (H, H q ) of two Jacobi opera- 
tors, where H is a short-range perturbation of a quasi-periodic finite-gap operator 
H q associated with the Riemann surface introduced in the previous section (see |20[ 
Ch. 9]). One key quantity is the transmission coefficient T{z). It is meromorphic 
in II + with finitely many simple poles in n+ n K precisely at the eigenvalues of the 
perturbed operator H. Since 

(3.1) |T(A)| 2 + |i?±(A)| 2 = l, AGE, 

it can be reconstructed from the reflection coefficients i?±(A) once we show how 
to reconstruct T(z) from its boundary values |T(A)| 2 = 1 - |i?±(A)| 2 , A e dU + . 
Rather than enter into more details here, see |lj (respectively [H]), we will focus 
on the reconstruction procedure only. 

We begin by deriving a formula for the Green function of 11+ : 

Lemma 3.1. The Green function o/II+ with pole at po is given by 

(3.2) g(z, z ) = -Re / uj PoPo , p = (z, +), p = (z , +), 

where po = po* (i.e., the complex conjugate on the other sheet) and oj pq is the 
normalized Abelian differential of the third kind with poles at p and q. 

Proof. First of all observe w PoPo = to Po E Q — ^PqEo an d set 

(3.3) uj poEo = r±(z, z )dz 

on II+. Since <-o Po e is continuous on the branch cuts, the corresponding values of 
r± must match up, that is, 

(3.4) limr+(A + is, zq) = limr_(A — is, zq), A e S. 

ej.0 ejO 

Moreover, 



(3.5) Wp oBo = r T (z, z )dz 

on Hence, 



(3.6) u poPa = lim(r+(A + ie, Zo) — r_(A — is, zo)d\ — 2iIm(r(A, zo))d\, A G S, 

ej.0 

where r(A, zq) = lim^o r +(^ + is, Zo), shows that oj p0tPo is purely imaginary on the 
boundary of 11+ . Together with the fact that the a-periods of to popo vanish this 
shows u popo is purely imaginary on <9II+. Hence g(z, z ) vanishes on dH + and 
since it has the proper singularity at zq by construction, we are done. □ 
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Clearly, we can extend g{z, Zq) to a holomorphic function on M\{zo} by dropping 
the real part. By abuse of notation we will denote this function by g{p,po) as well. 
However, note that g(p,po) will be multivalued with jumps in the imaginary part 
across ^-cycles. We will choose the path of integration in C\[i?o, E 2g+1 ] to guarantee 
a single- valued function. 

From the Green's function we obtain the Blaschke factor and the harmonic mea- 
sure (see e.g., |23j). Since we are mainly interested in the case where the poles are 
on the real line (since T(z) has all poles on the real line), we note the following 
relation which will be needed later on: 

Lemma 3.2. For p with ir{p) £ R\E we have 

rp rp 

(3.7) g(p,p)= u PP » = u PP *, 

J E J E(p) 

where E(p) is E if p < Eq, either E 2 j-i or E 2 j if p E (E 2 j-i, E 2 j), 1 < j < g, 
and E 2g+ i if p > E 2g+ i. 

Proof. By symmetry of the Green's function this holds at least when taking real 
parts. Since both quantities are real for ir(p) < E it holds everywhere. □ 

Now we come to the Blaschke factor 

pp 

(3.8) B(p, p) = exp (g(p, p) \ = exp ( / w pp , V n(p) € R, 

J E 

and first show that it can be written in terms of thcta functions. 
Lemma 3.3. The Blaschke factor is given by 

( , q] j3( n) e(A pa {p*) + a po {V)~^ po ) 9(A po (p)-A po (p)~a pa (V)+E pQ ) 
1 ' ' [P,P1 0(A po (p)+a po (V)-E po ) 9(A po (p)-A po (p*)-a po (V)+E po y 

where T> is any divisor of degree g—1 such thatT> p +T> andT> p * -\-T> are nonspecial. 
In addition, it satisfies 

(3.10) B{E ,p) = l, and B(p* , p) = B(p, p*) = B(p, p)" 1 ; 
it is real-valued for n(p) G (— oo, Eq). 

Proof. Both the Blaschke factor and the quotient of theta functions are multivalued 
meromorphic functions. Invoking the bilinear relations shows that we have the same 
jumps f~ Q around 6-cycles. Hence their quotient is single-valued. Moreover both 
have the same divisor T> p — T> p and hence the quotient is holomorphic and thus 
constant. 

To see that B{p* , p)B(j>, p) = 1, note that this function has no jumps and no 
poles and hence is constant. Since it is one at Eq it is one everywhere. □ 

Next we compute the harmonic measure of dll+. 

Lemma 3.4. The harmonic measure of dH+ with pole at A is given by 

(3.11) p(p, \)d\ = -Imw pB „(A) = --Im ( / u>\,o) d\, 

n 7T \J Ea J 

where ujx.o is the normalized Abelian differential of the second kind with a second- 
order pole at A. 
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Proof. All we have to do is to compute (2ir) 1 lim 6 |o Jj:S( 2 i ^ i 1£ ) (where the sign 
is chosen according to which side of £ one is interested) : 



(3.12) - |-Re 

ae 



X±ie 

= Imwpp = 2Imu pEo 



Eg 



e=0 



UJpE - 

The other formula follows similarly: 

21m / uj x ,o 



(3.13) - ^Re / 



6=0 



•lEo 

Note that 

-1 



= V^(A-^) +0(1) J dA 
for A near Eq and that the imaginary part has no singularity for A G <9I1 + . 
Now we can characterize the scattering data ([I], [5]): 

Theorem 3.5. Let T[z) be meromorphic in n + with simple poles at {pj}' =1 C R\E 
smc/i that T is continuous up to the boundary with the only possible simple zeros at 
the branch points. 

Then T(z) can be recovered from the boundary values ln|T(A)|, A £ S, via the 
Poisson- Jensen-type formula 

(3.14) T(z) = \ f[ B(p, pj)- 1 j exp (j- In \T\ 2 u pE ^j , p=(z, +), 
where we have identified pj with (pj,+) and defined T(p) = lim e |o T(X ± e), p = 

(A,±)ean+. 

Proof. The formula for T(z) holds by [25l Thm. 1], when taking absolute values. 
Since both sides are analytic and have equal absolute values, they can only differ 
by a constant of absolute value one. But both sides are positive at z — oo and 
hence this constant is one. □ 

Remark 3.6. A few remarks are in order: 

(i) The integrand in \3. 1J$ is not integrable at Eq and the integral has to be 
understood as a principal value. Otherwise, one can move the singularity 
away from <9II + which just alters the value by a constant. 

(ii) In scattering theory one has \T(p*)\ = \T(p)\, p S c?n +7 and under this 
assumption we have 



T(z) 




j exp^£ln(l- \R ± \ 2 )lu pp , 



where E(p) is defined in Lemma \3.2\ and the integral over £ is taken on 
the upper sheet. 
(iii) The Abelian differential is explicitly given by 

( Rljli + RljUp) <+ 2 + <+ 2 (g) , p , A 
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where P pq (z) is a polynomial of degree g — 1 which has to be determined 
from the normalization J a lo pp * — 0. In particular, 

In inverse scattering theory one uses (|3 . 1 [) to reconstruct T from the reflection 
coefficient R+ (or R_). Since neither the Blaschke factors nor the outer function 
in (|3.14p are in general single- valued on II + , we are naturally interested in when T 
is single- valued for given R±. 

Theorem 3.7. The transmission coefficient T defined via \3. 1J$ is single-valued if 
and only if the eigenvalues pj and the reflection coefficient R± satisfy 

(3.15) r T~- i ln(1 ~ |jR±|2)C£ € z - 

Proof. For T(z) to be single-valued we need lim e j T(x — ie) — lim £ | T(x + ie) 
for every i in a spectral gap. If x £ (-E^g-i, E21) is m the ^'th gap, the path of 
integration from Eq to A from above and back from A to Eq from below just gives 
the 6-cycle bg. Hence, 
(3.16) 

T(x + is) _ f A f 1 



S W^e) = ^ \h L "»* - M ' ln(1 " |i?±(A) ' } ' " A '° dA 



dn + J b e 



Evaluating the 6^-cycle using the usual bilinear relations finally yields 

(3.17) + =e*P| 2 ^E F'b- I H^-\R±?)Q 

eiO T{x -is) \ J^i Hi J 911+ 

and if the limit is supposed to be one, we are lead to (|3.15|) . □ 

The special case for an elliptic background with zero reflection coefficient was first 
obtained in [T3]. An analogous result was obtained in a different context by [2"T] . 
One should also emphasize that (|3.15[) is only a necessary condition for the solv- 
ability of the inverse scattering problem. Necessary and sufficient conditions are 
given in [4]. 

4. Trace formulas 

The transmission coefficient also plays a central role in the inverse scattering 
transform. Since it turns out to be the perturbation determinant of the pair (H, H q ), 
in the sense of Krein ([12j. |13|). its asymptotic expansion provides the conserved 
quantities of the Toda hierarchy ([5], [25], [23), 

■ 00 

(4.1) ]n r(z) = -X;-g T , r k = tr(H" - (H q )«). 



dz 

fe=i 



Relating this expansion with the one obtained by expanding (|3.14|) near z = 00, 
one obtains the usual trace formulas (also known as Case-type sum rules, [3]). 
Next, let LOj be the meromorphic differential 

(4.2) ^0 = <^oo + oo_ , Wfc = Woo+.-fc— 1 — Woo_;fc-1, 
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where u> p ,k is the Abelian differential of the second kind with a pole of order k + 2 
at p. Note that ujk is of the form 

(4.3) u k = ^Mdrr, 

where Pk(z) is a monic polynomial of degree g + k whose coefficients have to be 
determined from the fact that the a-cycles vanish and from the behavior at oo± 
(see [2D1 Eq. (13.30)]). 

Theorem 4.1. The following trace formulas are valid: 



ln(T(oc)) = - V ["' Woo+oo. + 4 / hi|T| 



LO 



0C-|_OO_ > 



(4.4) -±-r k+1 = Wfc + l/ln|T|^, 



where E{p) is defined in Lemma \3. 2\ and the integral over £ is taken on the upper 
sheet. 

Proof. By -jp;LOp( z )E |z=o = fc!^p ,fe-ii where z is a coordinate centered at poi we 
have 

(4-5) LUpEo = Woo + Bo +X^ 2fcw P' fc - 1 ' P =(->+)) 

and the claim follows. □ 



Acknowledgments. I want to thank Peter Yuditskii for several helpful discussions 
and for pointing out [2lJ [22J, [25] to me. In addition, I am also grateful to Fritz 
Gesztesy for discussions on this topic. 

References 

[1] N. I. Akhiezer, Elements of the Theory of Elliptic Functions, Amer. Math. Soc, Provi- 
dence, 1990. 

[2] W. Bulla, F. Gesztesy, H. Holden, and G. Teschl Algebro- Geometric Quasi-Periodic 
Finite-Gap Solutions of the Toda and Kac-van Moerbeke Hierarchies, Memoirs of the 
Amer. Math. Soc. 135/641, (1998). 

[3] K. M. Case, Orthogonal polynomials II, J. Math. Phys. 16, 1435-1440 (1975). 

[4] I. Egorova, J. Michor, and G. Teschl, Scattering Theory for Jacobi Operators with Quasi- 
Periodic Background, Comm. Math. Phys. 264-3, 811-842 (2006). 

[5] I. Egorova, J. Michor, and G. Teschl, Inverse scattering transform for the Toda hierarchy 
with guasi-periodic background, Proc. Amer. Math. Soc. (to appear). 

[6] H. Farkas and I. Kra, Riemann Surfaces, 2 n< ^ edition, GTM 71, Springer, New York, 
1992. 

[7] C. S. Gardner, J. M. Green, M. D. Kruskal, and R. M. Miura, A method for solving the 
Korteweg-de Vries eguation, Phys. Rev. Letters 19, 1095-1097 (1967). 

[8] F. Gesztesy and R. Svirsky, (m)KdV solitons on the background of guasi-periodic finite- 
gap solutions, Memoirs Amer. Math. Soc. 118, No. 563 (1995). 

[9] F. Gesztesy and G. Teschl, Commutation methods for Jacobi operators, J. Din. Eqs. 
128, 252-299 (1996). 

[10] F. Gesztesy, W. Schweiger, and B. Simon, Commutation methods applied to the mKdV- 

eguation, Trans. Amer. Math. Soc. 324, no. 2, 465-525 (1991). 
[11] R. Killip and B. Simon, Sum rules for Jacobi matrices and their applications to spectral 

theory, Ann. of Math. (2) 158, 253-321 (2003). 



s 



G. TESCHL 



[12] I. C. Gohberg and M. G. Krem, Introduction to the theory of linear nonselfadjoint 
operators, Transl. of Math. Mon. 18, Amer. Math. Soc, Providence, R.I. 1969. 

[13] M.G. Krein, Perturbation determinants and a formula for the traces of unitary and 
self-adjoint operators, Soviet. Math. Dokl. 3, 707-710 (1962). 

[14] E.A. Kuznetsov and A.V. Mikhailov, Stability of stationary waves in nonlinear weakly 
dispersive media, Soviet Phys. JETP 40, No.5, 855-859 (1975). 

[15] A. Laptev, S. Naboko, and O. Safronov, On new relations between spectral properties of 
Jacobi matrices and their coefficients, Comm. Math. Phys. 241, no. 1, 91—110 (2003). 

[16] J. Michor and G. Teschl, Trace formulas for Jacobi operators in connection with scatter- 
ing theory for quasi-periodic background, in Proc. of the conference Operator Theory and 
Applications in Mathematical Physics - OTAMP2004, J. Janas, P. Kurasov, A. Laptev, 
S. Naboko, and G. Stolz (eds.), Oper. Theory Adv. Appl., Birkhauser, Basel, (to appear). 

[17] F. Nazarov, F. Peherstorfer, A. Volberg, and P. Yuditskii, On generalized sum rules for 
Jacobi matrices, Int. Math. Res. Not. 2005 no. 3, 155-186 (2005). 

[18] B. Simon and A. Zlatos, Sum rules and the Szego condition for orthogonal polynomials 
on the real line, Comm. Math. Phys. 242, no. 3, 393-423 (2003). 

[19] G. Teschl, Inverse scattering transform for the Toda hierarchy, Math. Nach. 202, 163— 
171 (1999). 

[20] G. Teschl, Jacobi Operators and Completely Integrable Nonlinear Lattices, Math. Surv. 

and Mon. 72, Amer. Math. Soc, Providence, R.I. 2000. 
[21] Ju. Ja. Tomcuk, Orthogonal polynomials on a system of intervals of the real axis, 

Har'kov. Gos. Univ. Uc. Zap. Vol. XXIX, 93-128 (1963). (Russian) 
[22] Ju. Ja. Tomcuk, Orthogonal polynomials on a given system of arcs of the unit circle, 

Sov. Math. Dokl. 4, 931-934, 1963. 
[23] M. Tsuji, Potential Theory in modern Functional Analysis, Maruzen, Tokyo, 1959. 
[24] A. Volberg and P. Yuditskii, On the inverse scattering problem for Jacobi Matrices with 

the Spectrum on an Interval, a finite systems of intervals or a Cantor set of positive 

length, Commun. Math. Phys. 226, 567-605 (2002). 
[25] V. Voichick and L. Zalcman, Inner and outer functions on Riemann surfaces, Proc. 

Amer. Math. Soc. 16, 1200-1204 (1965). 
[26] A. Zlatos, Sum rules for Jacobi matrices and divergent Lieb-Thirring sums, J. Funct. 

Anal. 225, no. 2, 371-382 (2005). 

Faculty of Mathematics, Nordbergstrasse 15, 1090 Wien, Austria, and International 
Erwin Schrodinger Institute for Mathematical Physics, Boltzmanngasse 9, 1090 Wien, 
Austria 

E-mail address: Gerald. Teschl9univie . ac . at 
URL: http: //www. mat .univie . ac . at/~gerald/ 



